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Abstract 

We report on a novel scheme to perform efficient simulations of Scanning Tunneling Microscopy 
(STM) of molecules weakly bonded to surfaces. Calculations are based on a tight binding (TB) 
technique including self-consistency for the molecule to predict STM imaging and spectroscopy. 
To palliate the lack of self-consistency in the tunneling current calculation, we performed first 
principles density-functional calculations to extract the geometrical and electronic properties of 
the system. In this way, we can include, in the TB scheme, the effects of structural relaxation 
upon adsorption on the electronic structure of the molecule. This approach is applied to the study 
of regioregular poly(3-dodecylthiophene) (F3DDT) polymer chains adsorbed on highly oriented 
pyrolytic graphite (BOFG). Results of spectroscopic calculations are discussed and compared with 
recently obtained experimental data. 
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I. INTRODUCTION 



The Scanning Tunneling Microscope is a remarkable tool to probe objects and materials 
at the atomic scale and has then widely been used to characterize the adsorption of various 
molecules on metallic or semiconducting surfaces[? ????]. Moreover, STM can be used 
to understand the relations between the electronic properties of materials and their struc- 
ture at atomic scale. Thus, the analysis of STM images allows, i) to study the structural 
parameters of the system; ii) to study the nature of the probed electronic states through 
the analysis of the bias dependent STM images. Also, complementary, Scanning Tunneling 
Spectroscopy (STS) measurements allow to scan the local density of states of the system[? 
? ] . However, the correspondence between topographic or spectroscopic measurements, and 
the underlying electronic structure is not always obvious. The bonding nature between the 
molecule and the substrate, the charge transfer or the interaction between severals molecules 
are examples of parameters that can modify the electronic properties of the system with 
respect to its isolated parts. There is therefore a crucial need for efficient theoretical tools 
to tackle these issues. Most of the models developed so far are based either on a simplified 
perturbative Bardeen approach[? ] or on the Landauer formalismf? ]. These methods 
have been successfully developed on the basis of ab initio or semi-empirical techniques on 
individual molecules [? ? ] or small supramolecular systems [? ] but the main restriction, at 
the present time, remains the computing time needed for those calculations. The question 
of the suitability of such methods for describing large molecular systems can then be debated. 

In this paper, an alternative model combining both ab initio and semi-empirical ap- 
proaches for the simulation of STM experiments is presented. This method is particularly 
well suited for the case of physisorbed molecules for which the effects of the adsorption on 
the electronic structure of the molecule remain weak. Thus, the tunneling current can be 
calculated within the tight binding approximation with no need for self-consistency. How- 
ever, information about the geometrical structure and the electronic properties of the system 
remains essential and will first be extracted from first principles calculations. This method 
will be discussed and compared with STS spectra recently obtained for regioregular poly(3- 
dodecylthiophene) (P3DDT) chains adsorbed on highly oriented pyrolytic graphite (HOPG). 
Indeed, 7r-conjugated semi-conducting polymers have emerged as a new class of materials 
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which present both self-assembly and novel electronic features[? ? ]. The physics of poly(3- 
alkylthiophenes), which stand as a generic model owing to their remarkable semi-crystalline 
properties and resulting high carrier mobilities [? ] can then directly be addressed by STM. 
Therefore, a recent study of regioregular P3DDT adsorbed on HOPG has been devoted to 
determine, at the local scale, the precise relationship between structural organization and 
electronic properties. Details about these experiments will be published elsewhere. The aim 
of this paper is not to give a complete theoretical analysis of this experimental study but to 
present a generic method that can be applied to various systems. We will then focus on an 
application of our theoretical model to the study of STS spectra on single defectless P3DDT 
chains adsorbed on a graphene surface. 

II. DESCRIPTION OF THE METHODS 

A. Tight Binding Calculation of the Tunneling Current 

STM experiments are simulated by calculating the current between a tip, that is assumed 
to have a pyramidal apex[? ], and a graphite surface, on which is physisorbed a single chain 
of P3DDT, as a function of the applied bias V g to the STM tip. The calculation is based on 
a tight binding formulation of the elastic scattering theory considering that the Hamiltonian 
of the system can be written as H = Hq + V where Ho is the Hamiltonian of the three 
uncoupled regions (the STM tip, the molecule, and the substrate) and V their couplingf? ]. 
The current is then given by 

W) = X £ I0ir(*)|i)| 2 

iesjer 

x if( e j ~ Pr) ~ f( e i ~ Vs)} tf(ei - ej), (1) 

where \j) and \i) are the eigenstates of H in the tip (T) and in the substrate (S), respectively, 
6j and €i being their energies. Equation ([TJ assumes that the eigenstates of H incident from 
the surface and partially transmitted to the tip are occupied by electrons up to the Fermi 
level fig, and symmetrically those incident from the tip are filled up to fix (/ is the Fermi- 
Dirac distribution function and Ht — Us = e Vg)- As shown in Refs.? and ? , Eq.((l} is a 
generalization of the Landauer formula. The scattering operator T(e) is expressed at each 



3 



energy e in a Green's function formalism [? ] as 

T(e) = V + VG(e)V, 

where G is the Green's function of the coupled system obtained from the Green's function 
of the isolated regions 

G (e) = lim(e - H + i 7 ) 

7— *o 

via the Dyson's equation 

G=[J-G V]- 1 Go. 

A small imaginary part 7 is added in the Green's functions to avoid the divergence for 
energies corresponding to the poles of the Hamiltonian. It also have a physical meaning, 
and reflects the broadening of the electronic states due to the coupling between the different 
parts. Considering that the interaction of the molecule with the two electrodes is small, we 
take 7 = 0.05 eV. 

For the coupling of the molecule with the tip and the surface, the hopping terms V a p 
between two atomic orbitals \a) and \/3) follow an exponential decay with the inter-atomic 
distance d. The prescription of Ref.? is assumed so that 

V&do), (2) 

where do is the sum of the covalent radii of the two atoms and V^Jdo) is the hopping 
integral calculated from Harrison's rules [? ] at distance do. 

To compute the Green's functions needed to calculate the tunneling current, the electronic 
structure of the system has first to be determined. Considering that the molecule is only 
weakly interacting with the surface and the tip, as it is the case in physisorption, one can 
reasonably assume that its electronic structure will not be strongly affected by adsorption. 
On the other hand, electronic correlations inside the molecule are important and have to be 
taken into account. Hence, the electronic structure of the molecule is calculated using a self- 
consistent TB method described by Krzeminski et al[? ? ] and electron-electron interactions 
are treated by a method derived from the orthodox theory [? ]. In order to take into account 
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charging effects (i.e., the electrons injected by the electrodes are charging the molecule), the 
molecule has to be considered as an electronic quantum dot [? ], characterized by affinity 
and ionization levels, ef(V g ) and e^(V^), that amount to the injection of one electron or one 
hole, respectively. The relevant gap for STM simulations is then given by the quasi-particle 
gap (i.e., the gap obtained from a many body approach and that is different to the Kohn- 
Sham eigenvalues) defined as the difference between Electron Affinity (EA) and Ionization 
Potential (IP). Moreover, the injection of an extra charge inside the molecule will also lead 
to a response of the charges inside the metallic electrodes. This extra charge will then be 
screened by the graphite surface and the STM tip. This phenomenon, that tends to reduce 
the gap of the molecule, is taken into account by the image charge methodf? ]. Finally, the 
dependence of the molecular levels with the applied bias is described by 

ef\V g ) = ef\V 9 = 0) - V eV 9 (3) 

where rj is the same constant for all the electronic statesf? ]. The parameter rj is deter- 
mined by the geometry of the junction and describes the potential's drop between the two 
electrodes. n can take values between and 1, the 0.5 value corresponding to a symmetric 
interaction of the molecule with the two leads, and the lower and upper limits standing for 
a molecule whose coupling to one electrode is predominant [? ]. This parameter (which is 
the only tunable parameter of our model), as well as the equilibrium energetic state of the 
junction (V g = 0) will determine the electronic states involved in the tunneling process. 

B. Ab initio calculations 

Although the previously described TB scheme is sufficient to perform STM simulations, 
ab initio calculations can provide important information about the structural and electronic 
properties of the system. Indeed, we will use first principles calculations to extract the real 
adsorption configuration of the molecule on the surface and to determine the effects of the 
adsorption on the electronic structure of the molecule. We will show how these information 
can directly be integrated into the TB model. 

The structural and electronic properties of the adsorbed polymer on the graphite surface 
have been calculated by ab initio density-functional theory (DFT) simulations using a 



5 



plane wave basis set treating the electron-ion interaction within the projector augmented 
wave method (PAW) and the Generalized Gradient Approximation (GGA) for describing 
exchange-correlation effects as implemented in the VASP code.[? ] The relaxed structure 
has first been obtained by energy minimization of the periodic unit of the polymer chain 
(corresponding to two thiophene rings) on a graphene layer. The molecule has been 
positioned at 0.35 nm from the surface, which corresponds to typical distance for it-it 
interactions. These calculations were performed taking into account that the dodecyl chains 
of the polymer are epitaxed onto the zigzag pattern of the honey comb lattice of HOPG, as 
already experimentally evidenced. [? ? ] This imposes that the periodicity of the molecule 
matches the one of the graphene layer, and we will then consider the periodic unit of the 
system represented on FigG] for the rest of the ab initio calculations. Thus, the considered 
unit cell comprises a large enough graphene layer to avoid interaction between molecules 
of adjacent cells in the perpendicular direction of the polymer chain. Finally, to reduce 
the computional cost, the periodic unit of P3DDT will be considered without the dodecyl 
chains as they do not contribute to the tunneling current. Indeed, alkyl chains, made of 
a bonds between carbon atoms, are very insulating as the HOMO of the chain is several 
eV below the Fermi level of the surface. Thus, the dodecyl chains' properties will not been 
probed for voltages generally applied in STM experiments. 

We now need to extract several parameters, for the TB model, as far as the electronic 
structure of the molecule is concerned. We start by determining what are the effects of 
the adsorption on the electronic structure of the molecule. As already mentioned, in the 
case of a weak bonding with the surface, these effects should be small. In order to check 
this point, we look to dispersionless features in the band structure of the combined system. 
Indeed, since the system is periodic in one direction of space, it is possible to extract the 
band structure of both isolated surface and molecule and compare them to the one of the 
complete system. The three band structures, plotted along the T-X direction in the first 
Brillouin zone, are represented in Figj2j The results confirm that the adsorption of P3DDT 
on HOPG is relatively weak and that it does not modify the electronic structure of the 
polymer except for a shift of all the molecular states so that there is almost an alignment 
of the midgap of the polymer with the Fermi level of the graphite surface. 
This also suggests that no significant charge transfer effects from the surface to the polymer 
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FIG. 1: (Color online) Relaxed structure of the periodic system (polymer adsorbed on a graphene 
layer) obtained by imposing that the dodecyl chains of the polymer (not shown) match the zigzag 
pattern of the graphite lattice. In the perpendicular direction, the periodicity of the surface is 
taken large enough to avoid interactions between the 7r-conjugated part of molecules in neighboring 
cells. Carbon, sulphur and hydrogen atoms are represented using gray, yellow and green colors, 
respectively. 

layer should be expected. Finally, this justifies that the electronic structure of the three 
regions of the tunneling junction can be calculated separately for the calculation of the 
current. One just has to rigidly shift all the molecular levels to keep the same relative 
position of occupied and unoccupied states of the polymer with respect to the Fermi level 
of the surface (with a work function of graphene of = 4.5 eV below vacuum level). 

However, in the TB method developed to calculate the tunneling current, a finite sized 
system is required. Since the electronic gap of ^-conjugated polymers decreases and, at some 
point, saturates when increasing the number of periodic units[? ], it is possible to determine 
what are the number of thiophene rings needed to replace an infinite long polymer chain by 
a finite one with similar quasi-particle gap. This can be done by calculating the evolution of 
the gap of the polymer with the number of repeatition units. Hence, for a set a finite sized 
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FIG. 2: Calculated band structures for (a) the graphene surface, (b) the isolated polymer and (c) 
the complete polymer and graphite system. On the latter, the HOMO and LUMO states of the 
polymer are identified by dotted lines. The zero of energy in each panel is settled by the respective 
Fermi level of each system. The adsorption of the polymer on the surface only leads to a rigid shift 
of all molecular states. 

molecules (2,4,6 and 8 periodic units), the quasi-particle gap is calculated as the difference 
between Electron Affinity (EA) and Ionization Potential (IP) defined as 
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EA = E[-l] - E[0] 
IP = E[0] - E[+l] 



(4) 



where E[— 1], E[0] and -E[+l] are the total energies of the negatively charged, neutral and 
positively charged species, respectively. The evolution of the quasi-particle gap as a function 
of the inverse of the number of repeating unit is presented in FiglHl The results obtained 
for the four first polythiophenes have then been fitted to predict the evolution for a large 
number of rings. Fig£J] shows that the quasi-particle gap of these compounds saturates for 
a large number of thiophene groups. One finds that above 18 rings (octadecathienyl) the 
value of the gap does not evolve significantly. This molecule will therefore be taken as a 
good approximation to simulate the infinite polymer chain. The value of the HOMO-LUMO 
gap (at the Kohn-Sham level of DFT) of the infinite polymer is 1.49 eV (see Fig0b), that 
should be compared to the molecular gap (EA — IP = 2.33 eV for octadecathienyl) in order 
to get a glimpse of charging effects in the molecular system. In practice, taking into account 
the size of the considered octadecathienyl molecule, the graphite surface will be simulated 
by a single graphene sheet using 150 A; points in a super-cell made of 504 atoms and only 
one 7r orbital by carbon atom will be kept (V pp7T = —2.75 eV, for the TB parametrization). 

One should keep in mind that the band structures calculations for the extended system 
can only been performed for a neutral configuration. This means that we are only able, 
at this stage, to locate the HOMO (occupied states) and the LUMO (unoccupied states) 
of the molecule with respect to the Fermi level of the graphene surface, mimicking in a 
situation where no current is passing through the molecule. Notwithstanding, as already 
mentioned, the relevant gap for STM measurements is the so-called quasi-particles gap. 
Therefore, for those STM simulations the key quantity is to position all the molecular 
levels with respect to the Electron Affinity and Ionization Potential of the molecule. This 
study can be done easily by considering isolated polythiophenes and look at the evolution 
of the position of the HOMO and IP levels (with respect to the vacuum level) with the 
number of thiophene rings. The inset of FigOl shows that not only charging effects tend to 
reduce the molecular gap but that they also shift downwards all occupied states by a value 
of AE = —0.3 eV. Assuming that the same change in the IP occurs when the polymer is 
adsorbed on the surface, then all occupied states have to be shifted by the same value AE. 
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FIG. 3: Evolution of the calculated quasi-particles gap (EA-IP) of poly(thiophene)s as a function 
of the inverse of the number of rings (A) and the corresponding fitting curve. Inset: evolution of 
HOMO (□) and IP (o) levels as a function of the inverse of the number of rings and their fitting 
curve. Calculated values have been fitted by a function a\ + n ° 2 a3 , where 01,02,03 are fitting 
parameters and n the number of thiophene rings. 

A scissor operator is then applied to all unoccupied states so that the energy difference 
between occupied and unoccupied ones corresponds to the quasi-particle gap which value 
is 2.33 eV for the octadecathienyl molecule. All this energy level scheme is schematically 
illustrated in FigEJt. 

The subsequent relevant effect for the modeling of the STM current is the screening of 
the injected charge by the metallic surface and STM tip. This phenomenon, that tends 
to reduce the molecular gap, is taken into account by the image charge method. The 
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renormalised last occupied molecular level and the first unoccupied one are now referred to 
as IP and EA. The new relevant energetic situation is depicted in Figl^i. 

From a practical point of view, in order to incorporate all those effects in the TB calcu- 
lation of the tunneling current, we have to follow five steps :(i) to determine the electronic 
structure of the octadecathienyl molecule with the self-consistent TB method described in 
Ref.? and ? ; (ii) to rigidly shift all the molecular states so that the position of the 
midgap of the molecule with respect to the Fermi level of the graphite surface corresponds 
to the situation extracted from the band structure plotted on Fig|2t; (iii) to apply a shift 
of AE = —0.3 to all occupied states so that the last occupied one corresponds to the IP of 
the molecule; (iv) to apply a scissor operator to all unoccupied states so that the molecular 
gap corresponds to the quasi-particle gap (the position of the first unoccupied state now 
corresponds to the EA of the molecule); (v) finally, to equally reduce the gap by the value 
extracted from the image charge method described in Ref.? (in the present case, this 
corresponds to a reduction of the levels of 0.7 eV). 

C. Experiments 

P3DDT films with sub-monolayer coverage were prepared as described in Ref.? . STM 
experiments were performed under ultrahigh- vacuum conditions (base pressure below 5.10 -11 
mbars) using a VT Omicron system and mechanically prepared Ptlr tips. Measurements 
were recorded at room temperature in the low-current mode (LC-STM) using an adequate 
preamplifier which operation required bias voltage (Vg) to be applied to the tip and sam- 
ple to be grounded. Spectroscopic data were acquired in the current imaging tunneling 
spectroscopy mode (CITS) in which topographic images, in the constant current mode, are 
recorded simultaneously with I(V g ) curves taken on a 100 x 100 points grid, feedback loop 
being disabled during spectroscopic acquisition. Bias voltage was swept starting from the 
topographic regulation set point. Conductance spectra dI/dV g were recorded using a lock- 
in technique, with a topographic feedback loop gain adjusted to avoid oscillations in the 
z-images. 
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III. RESULTS AND DISCUSSION 



As a first application of our model, I(V g ) curve has been calculated for the bare graphene 
surface. In these calculations, the STM tip is positioned at a height giving a value of the 
tunneling current of the order of the experimental one (some pA, for a tip height of about 8 
A above the surface in the simulations). The comparison with experimental measurements is 
shown on Fig|3h with a good agreement, both curves being consistent with the semi-metallic 
nature of graphite. A completely different behavior is found for the polymer chain for which 
I{V g ) curves and conductance spectra exhibit an extended plateau (or conductance gap) 
with zero current. Calculations have been performed for severals values of the rj parameter 
between 0.5 and 0. We remind that this parameter is used to describe how the potential 
falls through the junction due to the coupling of the molecule with the substrate and the 
STM tip. The value rj = 0.5 indicates that the molecule is equally interacting with the two 
electrodes, whereas the case rj = reflects a predominant coupling with the surface. Values 
between 0.5 < rj < 1 should be excluded since a stronger interaction with the STM tip 
than with the substrate seems unreasonable. In Fig|Sj3, the calculated I(V g ) curves for three 
different values of rj (0, 0.33 and 0.5) are shown (since the position of the tip is kept constant 
for all the calculations, the three presented curves can directly be compared to each other). 
The shape of these curves is very similar and a zero conductance plateau is obtained for 
each cases. However, the width of this plateau is directly related to the value of rj, which 
suggests that the abrupt increase of the tunneling current is due to the contribution of the 
molecular levels. Indeed, the shift of these levels with the applied bias is related to the 
potential's drop in the junction (Eq|SJ). Hence, for a high enough negative bias, electrons 
will be injected from the STM tip to the surface via the EA (FigEb) and first unoccupied 
levels, as we stand in the situation where we have a quasi continuum of states. Similarly, 
for a sufficiently strong positive bias, electrons will be injected from the surface to the tip 
via the IP and last occupied levels of the molecule (Fig Eh); the values of these voltages 
depending on the value of the rj parameter. 

Another important feature of these I(V g ) curves is their asymmetry as a much higher 
current is obtained for positive biases than for negative ones. Additionally, the voltage 
for which a tunneling current suddenly increases is not the same at positive and negative 
polarity. Differently to the work of Terada et al. on P3HT on silicon[? where a 
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valence band-based conduction takes place whatever the bias voltage polarity, we have 
demonstrated that the conductance gap for P3DDT on HOPG is directly related to the 
bandgap of the polymer, as also evidenced for polydiacetylene on the same substrate [? ]. 
The asymmetry of the curves is then explained by the difference of voltages to be applied 
to inject electrons via the occupied or unoccupied molecular states (FigHJj). 

The comparison between the calculated spectra and experimental one is given on Fig|St. 
A very good agreement has been obtained for rj = 0.33 which is consistent with a stronger 
interaction of the molecule with the substrate due to the fact that, in the low-current STM 
mode, the tip-sample distance is relatively large. On the other hand, a value of 77 different 
from is also expected since the bonding of the molecule with the substrate is weak which 
supports the absence of charge transfer effects predicted by our ab initio calculations. As 
a matter of fact, within the band gap (FigEfc), the density of states takes no finite value, 
leading to a conductance gap of ca. 1.7 eV from experimental curves. In other words, there 
are no charge transfer effects from the substrate to the polymer for the creation of polaronic 
states inside the semi-conducting gap. 

Finally, one can mention that, for quasi one dimensional 7r-conjugated systems, 7r-band edge 
singularities might be expected. Such features have been reported from STS spectra on 
PDA[? ], but the electronic structure of P3DDT, as revealed by experimental and simulated 
curves, appears more complex. On measured spectra, bumps are observed on the falling edge 
of the last occupied molecular state instead of a well define single peak. dl/dV calculated 
at T = K shows, for its part, several peaks at positive voltages, in good agreement with 
the measured bumps if one takes into account thermal broadening effects. 

IV. CONCLUSION 

In conclusion, we have presented a generic method that combines first-principles and 
semi-empirical approaches for the simulation of STM experiments in the case of molecules 
weakly bonded to surfaces. A tight binding scheme has been developed for the calculation 
of tunneling current, whereas ab initio calculations have been used to extract information 
about the structural and electronic properties of the system. Although particularly adapted 
to long periodic molecules, such as polymer chains, it would be possible, in principle, to 
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apply the same methodology to smaller and non periodic ones. We have applied these 
methods to the study of defectless poly(3-dodecylthiophene) chains adsorbed on graphite, 
and showed that STS spectra can directly be related to the electronic structure of the 
polymer itself. Moreover, the comparison with experimental data has presented a very 
good agreement and supports the validity of our theoretical approach. Furthermore, the 
wide number of experimental results published by the scientific community on 7r-conjugated 
systems adsorbed on HOPG, highlights the importance of this theoretical work than can 
easily be applied to this great variety of systems. 
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FIG. 4: Energetic schemes of the tunneling junction, (a) Energy levels of octadecathienyl molecule 
with respect to the Fermi levels of the electrodes at zero bias : single-particle levels for the HOMO 
and LUMO, ionization levels corresponding to the charging of the molecule with one electron (EA) 
or one hole (IP) for the isolated molecule and for the molecule coupled to the electrodes ( EA 
and IP, respectively), (b), (c) Evolution of the bias-dependent molecular level positions taking 
i] = 0.33 (see text) : for a high enough negatjye bias, electrons are injected from the tip to the 
surface via EA level of the molecule, whereas for a sufficiently positive bias electrons are injected 
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I. INTRODUCTION 



The Scanning Tunneling Microscope is a remarkable tool to probe objects and materials 
at the atomic scale and has then widely been used to characterize the adsorption of various 
molecules on metallic or semiconducting surfaces [?????]. Moreover, STM can be used 
to understand the relations between the electronic properties of materials and their struc- 
ture at atomic scale. Thus, the analysis of STM images allows, i) to study the structural 
parameters of the system; ii) to study the nature of the probed electronic states through 
the analysis of the bias dependent STM images. Also, complementary, Scanning Tunneling 
Spectroscopy (STS) measurements allow to scan the local density of states of the system [? 
? ] . However, the correspondence between topographic or spectroscopic measurements, and 
the underlying electronic structure is not always obvious. The bonding nature between the 
molecule and the substrate, the charge transfer or the interaction between severals molecules 
are examples of parameters that can modify the electronic properties of the system with 
respect to its isolated parts. There is therefore a crucial need for efficient theoretical tools 
to tackle these issues. Most of the models developed so far are based either on a simplified 
perturbative Bardeen approach [? ] or on the Landauer formalism [? ]. These methods 
have been successfully developed on the basis of ab initio or semi-empirical techniques on 
individual molecules [? ? ] or small supramolecular systems [? ] but the main restriction, 
at the present time, remains the computing time needed for those calculations. The question 
of the suitability of such methods for describing large molecular systems can then be debated. 

In this paper, an alternative model combining both ab initio and semi-empirical ap- 
proaches for the simulation of STM experiments is presented. This method is particularly 
well suited for the case of physisorbed molecules for which the effects of the adsorption on 
the electronic structure of the molecule remain weak. Thus, the tunneling current can be 
calculated within the tight binding approximation with no need for self-consistency. How- 
ever, information about the geometrical structure and the electronic properties of the system 
remains essential and will first be extracted from first principles calculations. This method 
will be discussed and compared with STS spectra recently obtained for regioregular poly(3- 
dodecylthiophene) (P3DDT) chains adsorbed on highly oriented pyrolytic graphite (HOPG). 
Indeed, 7r-conjugated semiconducting polymers have emerged as a new class of materials 
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which present both self-assembly and novel electronic features [? ? ]. The physics of poly(3- 
alkylthiophenes), which stand as a generic model owing to their remarkable semi-crystalline 
properties and resulting high carrier mobilities [? ] can then directly be addressed by STM. 
Therefore, a recent study of regioregular P3DDT adsorbed on HOPG has been devoted to 
determine, at the local scale, the precise relationship between structural organization and 
electronic properties. Details about these experiments will be published elsewhere. The aim 
of this paper is not to give a complete theoretical analysis of this experimental study but to 
present a generic method that can be applied to various systems. We will then focus on an 
application of our theoretical model to the study of STS spectra on single defectless P3DDT 
chains adsorbed on a graphene surface. 

II. DESCRIPTION OF THE METHODS 

A. Tight Binding Calculation of the Tunneling Current 

STM experiments are simulated by calculating the current between a tip, that is assumed 
to have a pyramidal apex [? ], and a graphite surface, on which is physisorbed a single chain 
of P3DDT, as a function of the applied bias V g to the STM tip. The calculation is based on 
a tight binding formulation of the elastic scattering theory considering that the Hamiltonian 
of the system can be written as H = Hq + V where Hq is the Hamiltonian of the three 
uncoupled regions (the STM tip, the molecule, and the substrate) and V their coupling [? 
]. The current is then given by 

W) = X £ I0ir(*)|i)| 2 

iesjer 

x if( e j ~ Pr) ~ f( e i ~ Vs)} tf(ei - ej), (1) 

where \j) and \i) are the eigenstates of H in the tip (T) and in the substrate (S), respectively, 
6j and €i being their energies. Equation ([TJ assumes that the eigenstates of H incident from 
the surface and partially transmitted to the tip are occupied by electrons up to the Fermi 
level fig, and symmetrically those incident from the tip are filled up to fix (/ is the Fermi- 
Dirac distribution function and \xt — (J>s — e Vg)- As shown in Refs. ? and ? , Eq.(|l]) is a 
generalization of the Landauer formula. The scattering operator T(e) is expressed at each 
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energy e in a Green's function formalism [? ] as 

T(e) = V + VG(e)V, 

where G is the Green's function of the coupled system obtained from the Green's function 
of the isolated regions 

G (e) = lim(e - H + i 7 ) 

7— *o 

via the Dyson's equation 

G=[J-G V]- 1 Go. 

A small imaginary part 7 is added in the Green's functions to avoid the divergence for 
energies corresponding to the poles of the Hamiltonian. It also have a physical meaning, 
and reflects the broadening of the electronic states due to the coupling between the different 
parts. Considering that the interaction of the molecule with the two electrodes is small, we 
take 7 = 0.05 eV. 

For the coupling of the molecule with the tip and the surface, the hopping terms V a p 
between two atomic orbitals \a) and \j3) follow an exponential decay with the inter-atomic 
distance d. The prescription of Ref. ? is assumed so that 

V&do), (2) 

where do is the sum of the covalent radii of the two atoms and V^Jdo) is the hopping 
integral calculated from Harrison's rules [? ] at distance do. 

To compute the Green's functions needed to calculate the tunneling current, the electronic 
structure of the system has first to be determined. Considering that the molecule is only 
weakly interacting with the surface and the tip, as it is the case in physisorption, one can 
reasonably assume that its electronic structure will not be strongly affected by adsorption. 
On the other hand, electronic correlations inside the molecule are important and have to be 
taken into account. Hence, the electronic structure of the molecule is calculated using a self- 
consistent TB method described by Krzeminski et al [? ? ] and electron-electron interactions 
are treated by a method derived from the orthodox theory [? ] . In order to take into account 
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Vapid) = exp 



d 



charging effects (i.e., the electrons injected by the electrodes are charging the molecule), the 
molecule has to be considered as an electronic quantum dot [? ], characterized by affinity 
and ionization levels, ef(V g ) and e^(V^), that amount to the injection of one electron or one 
hole, respectively. The relevant gap for STM simulations is then given by the quasi-particle 
gap (i.e., the gap obtained from a many body approach and that is different to the Kohn- 
Sham eigenvalues) defined as the difference between Electron Affinity (EA) and Ionization 
Potential (IP). Moreover, the injection of an extra charge inside the molecule will also lead 
to a response of the charges inside the metallic electrodes. This extra charge will then be 
screened by the graphite surface and the STM tip. This phenomenon, that tends to reduce 
the gap of the molecule, is taken into account by the image charge method [? ]. Finally, the 
dependence of the molecular levels with the applied bias is described by 

ef\V g ) = ef\V 9 = 0) - V eV 9 (3) 

where 77 is the same constant for all the electronic states [? ]. The parameter rj is deter- 
mined by the geometry of the junction and describes the potential's drop between the two 
electrodes. 77 can take values between and 1, the 0.5 value corresponding to a symmetric 
interaction of the molecule with the two leads, and the lower and upper limits standing for 
a molecule whose coupling to one electrode is predominant [? ]. This parameter (which is 
the only tunable parameter of our model), as well as the equilibrium energetic state of the 
junction (V g = 0) will determine the electronic states involved in the tunneling process. 

B. Ab initio calculations 

Although the previously described TB scheme is sufficient to perform STM simulations, 
ab initio calculations can provide important information about the structural and electronic 
properties of the system. Indeed, we will use first principles calculations to extract the real 
adsorption configuration of the molecule on the surface and to determine the effects of the 
adsorption on the electronic structure of the molecule. We will show how these information 
can directly be integrated into the TB model. 

The structural and electronic properties of the adsorbed polymer on the graphite surface 
have been calculated by ab initio density-functional theory (DFT) simulations using a 
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plane wave basis set treating the electron-ion interaction within the projector augmented 
wave method (PAW) and the Generalized Gradient Approximation (GGA) for describing 
exchange-correlation effects as implemented in the VASP code [? ]. The relaxed structure 
has first been obtained by energy minimization of the periodic unit of the polymer chain 
(corresponding to two thiophene rings) on a graphene layer. The molecule has been 
positioned at 0.35 nm from the surface, which corresponds to typical distance for ir-ir 
interactions. These calculations were performed taking into account that the dodecyl chains 
of the polymer are epitaxed onto the zigzag pattern of the honey comb lattice of HOPG, as 
already experimentally evidenced [? ? ]. This imposes that the periodicity of the molecule 
matches the one of the graphene layer, and we will then consider the periodic unit of the 
system represented on Fig.l for the rest of the ab initio calculations. Thus, the considered 
unit cell comprises a large enough graphene layer to avoid interaction between molecules 
of adjacent cells in the perpendicular direction of the polymer chain. Finally, to reduce 
the computional cost, the periodic unit of P3DDT will be considered without the dodecyl 
chains as they do not contribute to the tunneling current. Indeed, alkyl chains, made of 
a bonds between carbon atoms, are very insulating as the HOMO of the chain is several 
eV below the Fermi level of the surface. Thus, the dodecyl chains' properties will not been 
probed for voltages generally applied in STM experiments. 

We now need to extract several parameters, for the TB model, as far as the electronic 
structure of the molecule is concerned. We start by determining what are the effects of 
the adsorption on the electronic structure of the molecule. As already mentioned, in the 
case of a weak bonding with the surface, these effects should be small. In order to check 
this point, we look to dispersionless features in the band structure of the combined system. 
Indeed, since the system is periodic in one direction of space, it is possible to extract the 
band structure of both isolated surface and molecule and compare them to the one of the 
complete system. The three band structures, plotted along the T-X direction in the first 
Brillouin zone, are represented in Fig. 2. The results confirm that the adsorption of P3DDT 
on HOPG is relatively weak and that it does not modify the electronic structure of the 
polymer except for a shift of all the molecular states so that there is almost an alignment 
of the midgap of the polymer with the Fermi level of the graphite surface. 
This also suggests that no significant charge transfer effects from the surface to the polymer 



6 



layer should be expected. Finally, this justifies that the electronic structure of the three 
regions of the tunneling junction can be calculated separately for the calculation of the 
current. One just has to rigidly shift all the molecular levels to keep the same relative 
position of occupied and unoccupied states of the polymer with respect to the Fermi level 
of the surface (with a work function of graphene of = 4.5 eV below vacuum level). 

However, in the TB method developed to calculate the tunneling current, a finite sized 
system is required. Since the electronic gap of ^-conjugated polymers decreases and, at some 
point, saturates when increasing the number of periodic units [? ], it is possible to determine 
what are the number of thiophene rings needed to replace an infinite long polymer chain by 
a finite one with similar quasi-particle gap. This can be done by calculating the evolution 
of the gap of the polymer with the number of repetition units. Hence, for a set a finite sized 
molecules (2,4,6 and 8 periodic units), the quasi-particle gap is calculated as the difference 
between Electron Affinity (EA) and Ionization Potential (IP) defined as 

EA = E[-1] -E[0] 

IP = E[0] - E[+l] (4) 

where E[— 1], E[0] and are the total energies of the negatively charged, neutral and 

positively charged species, respectively. The evolution of the quasi-particle gap as a function 
of the inverse of the number of repeating unit is presented in Fig. 3. The results obtained 
for the four first polythiophenes have then been fitted to predict the evolution for a large 
number of rings. Fig. 3 shows that the quasi-particle gap of these compounds saturates for 
a large number of thiophene groups. One finds that above 18 rings (octadecathienyl) the 
value of the gap does not evolve significantly. This molecule will therefore be taken as a 
good approximation to simulate the infinite polymer chain. The value of the HOMO-LUMO 
gap (at the Kohn-Sham level of DFT) of the infinite polymer is 1.49 eV (see Fig. 2 b), that 
should be compared to the molecular gap (EA — IP = 2.33 eV for octadecathienyl) in order 
to get a glimpse of charging effects in the molecular system. In practice, taking into account 
the size of the considered octadecathienyl molecule, the graphite surface will be simulated 
by a single graphene sheet using 150 A; points in a super-cell made of 504 atoms and only 
one 7r orbital by carbon atom will be kept (V pp7T = —2.75 eV, for the TB parametrization) . 
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One should keep in mind that the band structures calculations for the extended system 
can only been performed for a neutral configuration. This means that we are only able, 
at this stage, to locate the HOMO (occupied states) and the LUMO (unoccupied states) 
of the molecule with respect to the Fermi level of the graphene surface, mimicking in a 
situation where no current is passing through the molecule. Notwithstanding, as already 
mentioned, the relevant gap for STM measurements is the so-called quasi-particles gap. 
Therefore, for those STM simulations the key quantity is to position all the molecular 
levels with respect to the Electron Affinity and Ionization Potential of the molecule. This 
study can be done easily by considering isolated polythiophenes and look at the evolution 
of the position of the HOMO and IP levels (with respect to the vacuum level) with the 
number of thiophene rings. The inset of Fig. 3 shows that not only charging effects tend to 
reduce the molecular gap but that they also shift downwards all occupied states by a value 
of AE = —0.3 eV. Assuming that the same change in the IP occurs when the polymer is 
adsorbed on the surface, then all occupied states have to be shifted by the same value AE. 
A scissor operator is then applied to all unoccupied states so that the energy difference 
between occupied and unoccupied ones corresponds to the quasi-particle gap which value 
is 2.33 eV for the octadecathienyl molecule. All this energy level scheme is schematically 
illustrated in Fig. 4 a. 

The subsequent relevant effect for the modeling of the STM current is the screening of 
the injected charge by the metallic surface and STM tip. This phenomenon, that tends 
to reduce the molecular gap, is taken into account by the image charge method. The 
renormalised last occupied molecular level and the first unoccupied one are now referred to 
as IP and EA. The new relevant energetic situation is depicted in Fig. 4 a. 

From a practical point of view, in order to incorporate all those effects in the TB calcu- 
lation of the tunneling current, we have to follow five steps :(i) to determine the electronic 
structure of the octadecathienyl molecule with the self-consistent TB method described in 
Ref. ? and ? ; (ii) to rigidly shift all the molecular states so that the position of the 
midgap of the molecule with respect to the Fermi level of the graphite surface corresponds 
to the situation extracted from the band structure plotted on Fig. 2 c; (iii) to apply a shift 
of AE = —0.3 to all occupied states so that the last occupied one corresponds to the IP of 
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the molecule; (iv) to apply a scissor operator to all unoccupied states so that the molecular 
gap corresponds to the quasi-particle gap (the position of the first unoccupied state now 
corresponds to the EA of the molecule); (v) finally, to equally reduce the gap by the value 
extracted from the image charge method described in Ref. ? (in the present case, this 
corresponds to a reduction of the levels of 0.7 eV). 

C. Experiments 

P3DDT films with sub-monolayer coverage were prepared as described in Ref. ? . STM 
experiments were performed under ultrahigh- vacuum conditions (base pressure below 5.10~ n 
mbars) using a VT Omicron system and mechanically prepared Ptlr tips. Measurements 
were recorded at room temperature in the low-current mode (LC-STM) using an adequate 
preamplifier which operation required bias voltage (Vg) to be applied to the tip and sam- 
ple to be grounded. Spectroscopic data were acquired in the current imaging tunneling 
spectroscopy mode (CITS) in which topographic images, in the constant current mode, are 
recorded simultaneously with I(V g ) curves taken on a 100 x 100 points grid, feedback loop 
being disabled during spectroscopic acquisition. Bias voltage was swept starting from the 
topographic regulation set point. Conductance spectra dI/dV g were recorded using a lock- 
in technique, with a topographic feedback loop gain adjusted to avoid oscillations in the 
z-images. 

III. RESULTS AND DISCUSSION 

As a first application of our model, I(V g ) curve has been calculated for the bare graphene 
surface. In these calculations, the STM tip is positioned at a height giving a value of the 
tunneling current of the order of the experimental one (some pA, for a tip height of about 8 
A above the surface in the simulations). The comparison with experimental measurements is 
shown on Fig. 5 a with a good agreement, both curves being consistent with the semi-metallic 
nature of graphite. A completely different behavior is found for the polymer chain for which 
I(V g ) curves and conductance spectra exhibit an extended plateau (or conductance gap) 
with zero current. Calculations have been performed for severals values of the r] parameter 
between 0.5 and 0. We remind that this parameter is used to describe how the potential 
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falls through the junction due to the coupling of the molecule with the substrate and the 
STM tip. The value rj = 0.5 indicates that the molecule is equally interacting with the two 
electrodes, whereas the case rj = reflects a predominant coupling with the surface. Values 
between 0.5 < r\ < 1 should be excluded since a stronger interaction with the STM tip than 
with the substrate seems unreasonable. In Fig. 5 b, the calculated I(V g ) curves for three 
different values of rj (0, 0.33 and 0.5) are shown (since the position of the tip is kept constant 
for all the calculations, the three presented curves can directly be compared to each other). 
The shape of these curves is very similar and a zero conductance plateau is obtained for each 
cases. However, the width of this plateau is directly related to the value of rj, which suggests 
that the abrupt increase of the tunneling current is due to the contribution of the molecular 
levels. Indeed, the shift of these levels with the applied bias is related to the potential's drop 
in the junction (EqJ3J). Hence, for a high enough negative bias, electrons will be injected 
from the STM tip to the surface via the EA (Fig. 4 b) and first unoccupied levels, as we 
stand in the situation where we have a quasi continuum of states. Similarly, for a sufficiently 
strong positive bias, electrons will be injected from the surface to the tip via the IP and 
last occupied levels of the molecule (Fig. 4 c); the values of these voltages depending on the 
value of the rj parameter. 

Another important feature of these I(V g ) curves is their asymmetry as a much higher 
current is obtained for positive biases than for negative ones. Additionally, the voltage 
for which a tunneling current suddenly increases is not the same at positive and negative 
polarity. Differently to the work of Terada et al. on P3HT on silicon [? where a 
valence band-based conduction takes place whatever the bias voltage polarity, we have 
demonstrated that the conductance gap for P3DDT on HOPG is directly related to the 
bandgap of the polymer, as also evidenced for polydiacetylene on the same substrate [? ]. 
The asymmetry of the curves is then explained by the difference of voltages to be applied 
to inject electrons via the occupied or unoccupied molecular states (Fig. 4 c). 

The comparison between the calculated spectra and experimental one is given on Fig. 5 c. 
A very good agreement has been obtained for rj = 0.33 which is consistent with a stronger 
interaction of the molecule with the substrate due to the fact that, in the low-current STM 
mode, the tip-sample distance is relatively large. On the other hand, a value of rj different 
from is also expected since the bonding of the molecule with the substrate is weak which 
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supports the absence of charge transfer effects predicted by our ab initio calculations. As 
a matter of fact, within the band gap (Fig. 5 c), the density of states takes no finite value, 
leading to a conductance gap of ca. 1.7 eV from experimental curves. In other words, there 
are no charge transfer effects from the substrate to the polymer for the creation of polaronic 
states inside the semi-conducting gap. 

Finally, one can mention that, for quasi one dimensional 7r-conjugated systems, 7r-band edge 
singularities might be expected. Such features have been reported from STS spectra on 
PDA[? ], but the electronic structure of P3DDT, as revealed by experimental and simulated 
curves, appears more complex. On measured spectra, bumps are observed on the falling edge 
of the last occupied molecular state instead of a well define single peak. dl/dV calculated 
at T = K shows, for its part, several peaks at positive voltages, in good agreement with 
the measured bumps if one takes into account thermal broadening effects. 

IV. CONCLUSION 

In conclusion, we have presented a generic method that combines first-principles and 
semi-empirical approaches for the simulation of STM experiments in the case of molecules 
weakly bonded to surfaces. A tight binding scheme has been developed for the calculation 
of tunneling current, whereas ab initio calculations have been used to extract information 
about the structural and electronic properties of the system. Although particularly adapted 
to long periodic molecules, such as polymer chains, it would be possible, in principle, to apply 
the same methodology to smaller and non periodic ones. We have applied these methods to 
the study of defectless poly(3-dodecylthiophene) chains adsorbed on graphite, and showed 
that STS spectra can directly be related to the electronic structure of the polymer itself. 
The comparison with experimental data has presented a very good agreement and supports 
the validity of our theoretical approach. Furthermore, the wide number of experimental 
results published by the scientific community on 7r-conjugated systems adsorbed on HOPG, 
highlights the importance of this theoretical work than can easily be applied to this great 
variety of systems. 
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Figures caption : 

Figure 1. (Color online) Relaxed structure of the periodic system (polymer adsorbed 
on a graphene layer) obtained by imposing that the dodecyl chains of the polymer (not 
shown) match the zigzag pattern of the graphite lattice. In the perpendicular direction, 
the periodicity of the surface is taken large enough to avoid interactions between the tt- 
conjugated part of molecules in neighboring cells. Carbon, sulfur and hydrogen atoms are 
represented using gray, yellow and green colors, respectively. 

Figure 2. Calculated band structures for (a) the graphene surface, (b) the isolated 
polymer and (c) the complete polymer and graphite system. On the latter, the HOMO and 
LUMO states of the polymer are identified by dotted lines. The zero of energy in each panel 
is settled by the respective Fermi level of each system. The adsorption of the polymer on 
the surface only leads to a rigid shift of all molecular states. 

Figure 3. Evolution of the calculated quasi-particles gap (EA-IP) of poly(thiophene)s 
as a function of the inverse of the number of rings (A) and the corresponding fitting curve. 
Inset: evolution of HOMO (□) and IP (o) levels as a function of the inverse of the number 
of rings and their fitting curve. Calculated values have been fitted by a function a% + 
where a±, a,2, are fitting parameters and n the number of thiophene rings. 

Figure 4. Energetic schemes of the tunneling junction, (a) Energy levels of octade- 
cathienyl molecule with respect to the Fermi levels of the electrodes at zero bias : single- 
particle levels for the HOMO and LUMO, ionization levels corresponding to the charging 
of the molecule with one electron (EA) or one hole (IP) for the isolated molecule and for 
the molecule coupled to the electrodes ( EA and IP, respectively), (b), (c) Evolution of 
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the bias-dependent molecular level positions taking 77 = 0.33 (see text) : for a high enough 
negative bias, electrons are injected from the tip to the surface via EA level of the molecule, 
whereas for a sufficiently positive bias electrons are injected form the surface to the tip via 
the IP level of the molecule. 

Figure 5. (Color Online) (a) Calculated (plain lines) and experimental (dotted lines) 
I(V g ) curves for the bare substrate (red) and the polymer chain (black). The calculated 
curve for the octadecathienyl molecule is presented for r\ = 0.33 (see text), (b) Calculated 
I(V g ) curves of the polymer chain for different values of rj. (c) Theoretical and experimental 
conductance dI/dV g spectra obtained for P3DDT. The presented calculated spectra has been 
obtained for 77 = 0.33. 
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M. Dubois et al., Figure 2 
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M. Dubois et al., Figure 3 
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M. Dubois et al., Figure 5 
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